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Abstract 

We present measurements of the branching fraction and time-dependent CP asymmetries in 

J/iIj TfO decays based on (231.8 ± 2.6) x 10'^ r(4S)^ BB decays collected with the BABAR 
detector at the PEP-II asymmetric-energy B factory at SLAC during the years 1999-2004. We 
obtain a branching fraction B{B^ J/i/j-k^) = (1.94 ± 0.22 (stat) ± 0.17 (syst))x 10"^ We also 
measure the CP asymmetry parameters C = —0.21 it 0.26 (stat) it 0.09 (syst) and S = —0.68 it 
0.30 (stat) lb 0.04 (syst). All results presented in this paper are preliminary. 
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1 Introduction 



Charge conjugation-parity (CP) violation in the B meson system has been estabhshed by the 
BABAR [1] and Belle [2] collaborations. The Standard Model (SM) of electroweak interactions 
describes CP violation as a consequence of a complex phase in the three-generation Cabibbo- 
Kobayashi-Maskawa (CKM) quark-mixing matrix [3]. Measurements of CP asymmetries in the 
proper-time distribution of neutral B decays to CP eigenstates containing a charmonium and 
meson provide a precise measurement of sin 2/3 [4], where (3 is arg [ — V^d^cb/^d^b ] ^^'^ ^ij 
CKM matrix elements. 

The decay B^ — ^ J/i/jtt^ is a CP-even Cabibbo-suppressed b — ccd transition for which, in the 
absence of loop (penguin) amplitudes, the SM predicts that the two CP asymmetry coefficients, 
S, the interference between mixing and decay, and C, the direct CP asymmetry, are S = — sin 2/3, 
and C = 0. S and C are defined as: 

22mA ^_ 1 - I^P 

__ and ^=TTW 



S^TTJu^ and C^^^, (1) 



where A is a complex parameter that depends on both the B^-B^ oscillation amplitude and the 
amplitudes describing B^ and B^ decays to the J/ipir^ final state. The tree and penguin amplitudes 
expected to dominate this decay are shown in Figure 1. 




Figure 1: Feynman diagrams of the color suppressed tree (left) and gluonic penguin (right) ampli- 
tudes contributing to the B'^ —>■ J/ijj 7r° decay. 



The b ^ ccd tree amplitude has the same weak phase as the b ^ ccs modes (e.g. the CP-odd 
decay P° J/i^Kl). The b ^ ccd penguin amplitudes have a different weak phase than the tree 
amplitude. If there is a significant penguin amplitude in B^ — > J/ij) t:^ , then one will measure a 
value of S that differs from — sin 2/3, and a value of C that differs from zero [5] . 

In this paper, we present an update of previous BABMi measurements of the branching fraction 
and time-dependent CP violating {CPV) asymmetries in B^ J if' [6, 7]. The preliminary 
results presented here are obtained using 210.6 fb~^ of data. BABPR and Belle have both previously 
presented measurements of the B^ ^ J/il) 7r° branching fraction using T(4S)— >^ BB decays. These 
are: 
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BABAR: (2.0 ± 0.6 (stat) ± 0.2 (syst))xlO-s (from 20.7 fb^i) [6], 
Belle: (2.3 ± 0.5 (stat) ± 0.2 (syst))xlO-^ (from 29.4 fb-^) [8]. 



Both the BABAR and Belle notations in denoting the magnitude of the direct CP asymmetry, 
where C {BABAR) = —A (Belle), the previous measurements of S and C (A) are: 

BABAR: S = 0.05 ± 0.49 ± 0.16 C = 0.38 ± 0.41 ± 0.09 (from 81.1 fb^^) [7], 
Belle: S = -0.72 ± 0.42 ± 0.09 A = -0.01 ± 0.29 ± 0.03 (from 140.0 fb^^) [9]. 



2 The BaBAR detector and dataset 

The data used in this analysis were collected with the BABAR detector at the PEP-II asymmetric 
e'''e~storage ring from 1999 to 2004. This represents a total integrated luminosity of 210.6 fb~^ 
taken at the T(4S) resonance (onpeak), corresponding to a sample of 231.8 ± 2.6 million BB pairs. 
An additional 21.6 fb~^ of data, collected at approximately 40 MeV below the T(4S) resonance, is 
used to study background from e'^e~ qq {q = u, d, s, c) continuum events. 

The BABAR detector is described elsewhere [10]. Surrounding the interaction point is a 5 layer 
double-sided silicon vertex tracker (SVT) which provides precise reconstruction of track angles and 
B decay vertices. A 40 layer drift chamber (DCH) surrounds the SVT and provides measurements 
of the transverse momenta for charged particles. Both the SVT and the DCH operate in a 1.5 
T solenoidal magnetic field. Charged hadron identification is achieved through measurements of 
particle energy-loss (dE/dx) in the tracking system and Cerenkov angle obtained from a detector 
of internally reflected Cerenkov light (DIRC). This is surrounded by a segmented CsI(Tl) electro- 
magnetic calorimeter (EMC) which provides photon detection, electron identification, and is used 
to reconstruct neutral hadrons. Finally, the instrumented flux return (IFR) of the magnet allows 
discrimination of muons from pions. 



3 Analysis Method 

We study B^ J/^ vr'' decays in BB candidate events from combinations of J/tp £~^i~ (£ = 
e, fi) and tt" —>■ 77 candidates. A detailed description of the J/tp — >■ £~^£~ selection can be found 
elsewhere [6]. For the J/ip — > e+e" {J/tp — > jU+jU") channel, the invariant mass of the lepton pair 
is required to be between 3.06 and 3.12GeV/c2 (3.07 and 3.13GeV/c2). 

We form tt^ 77 candidates with an invariant mass 100 < m^^ < 160 MeV/c^ from pairs of 
photon candidates which have been identified as clusters in the EMC. These clusters are required to 
be isolated from any charged tracks, carry a minimum energy of 30 MeV, and have a lateral energy 
distribution consistent with that of a photon. Each 7r° candidate is required to have a minimum 
energy of 200 MeV and is constrained to the nominal mass [11]. 

The B^ J/tJj TT^ candidates (Brec) arc reconstructed from these J/tjj and tt^ candidates and 
constrained to originate from the e^e^ interaction point using a geometric fit. Finally, the J/ip 
and 7r° candidates are combined using 4-momentum addition. 

We use two kinematic variables, mES AE, in order to isolate the signal. 
TIES = \/(\/^/2)^ — (Pb)^ is the beam-energy substituted mass, where ^/s is the center-of-mass 
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(CM) energy, and therefore \fsl1 is the beam energy in the CM frame, is the -B-candidate 
momentum in the CM frame. IS.E is the difference between the S-candidate energy and the beam 
energy in the e+e" CM frame. We require mES > 5.2GeV/c^ and |A£'| < 0.3 GeV. 

A significant source of background is due to e^e~ qq (q = u,d,s,c) continuum events. 
We combine several kinematic and topological variables into a Fisher discriminant (F) to provide 
additional separation between signal and continuum. The three variables Lq, L2 and cos{6h) 
are inputs to F. Lq and L2 are the zeroth- and second-order Legendre polynomial moments; 
Lq = J2i I Pi I L2 = X^i I Pi I where p* are the CM momenta for the tracks and neutral 

calorimeter clusters that are not associated with the signal candidate. The 6i are the angles between 
p* and the thrust axis of the signal candidate. Oh is the angle between one of the leptons and the 
B candidate in the J/^ rest frame. 

We use multivariate algorithms to identify signatures of B decays that determine (tag) the flavor 
of the decay of the other B in the event (Btag) to be either a B^ or B^. The flavor tagging algorithm 
used is described in more detail elsewhere [12]. In brief, we define seven mutually exclusive tagging 
categories. These are (in order of decreasing signal purity) Lepton, Kaon 1, Kaon 2, Kaon-Pion, 
Pion, Other, and No-Tag. The total effective tagging efficiency of this algorithm is (30.5 ± 0.4)%. 

The decay rate /+ (/_) of 5° decays to a CP eigenstate, when Btag is a B^ (-8°), is: 

-\At\/TgO 

f±iAt) = — [1 ± 5sin(AmdAt) ^ Ccos(AmdAt)], (2) 

4tbo 

where At is the difference between the proper decay times of the Brec and Bfag mesons, r^o = 
1.536 ± 0.014 ps is the B° lifetime and Ama = (0.502 ± 0.007) xlO'^^s is the B^-B^ oscillation 
frequency [11]. The decay width difference between the B^ mass eigenstates is assumed to be zero. 

The time interval At is calculated from the measured separation Az between the decay vertices 
of Brec and Btag along the collision axis (z). The vertex of Brec is from the lepton tracks that come 
from the J/^l^ and the vertex of Btag is constructed from the remaining tracks in the event that do 
not belong to Brec, whilst using constraints from the beam spot location and the Brec momentum. 
We accept events with \At\ < 20 ps whose uncertainty is less than 2.5 ps. 

After all of the selection cuts mentioned above have been applied, the average multiplicity 
is approximately 1.1, indicating some events still have multiple candidates. In these events, we 
randomly choose one candidate to be used in the fit. After this step, the signal efficiency is 22.0% 
and a total of 1318 onpeak events are selected. 

In addition to signal and continuum background events, there are also B backgrounds present 
in the data after applying the selection cuts above. We divide the B backgrounds into the following 
types: (i) B^ ^ J/ij; Kg, (ii) generic neutral B meson decays, and (iii) generic charged B meson 
decays. Prom Monte Carlo (MC) we expect 153 ±9, 68 ±14 and 314 ±63 events of these background 
types, respectively. The generic neutral B meson decays do not include signal or B^ J/il) Kg 
events. The generic B decay backgrounds are dominated by contributions from B J /'ijjX (inclu- 
sive charmonium final states). In particular the generic charged B meson decay backgrounds are 
dominated by B^ — > J/i)p'^ decays. The B'^ ^ J/ip Kg background was studied separately since 
it is a well understood decay with respect to time-dependent analysis. 

Wc perform an extended unbinned maximum likelihood fit to the B candidate sample, where 
the discriminating variables used in the fit are mES) AE, F and At. The values of the signal yield, 
S and C are simultaneously extracted. 

The signal tties distribution is described by a Gaussian with an exponential tail [13]. We pa- 
rameterise the mES distribution for continuum and neutral generic B background with a phase 
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space distribution [14]. As there are significant correlations between mss and AE for the charged 
generic B background, we parameterise these variables with two-dimensional non-parametric prob- 
ability density functions (PDFs). We use two-dimensional non-parametric PDFs when describing 
the TJiES-AE distribution for J/tpKs. The AE distribution for signal events is modeled with 

a Gaussian with an exponential tail on the negative side to model energy leakage in the EMC, plus 
a polynomial contribution. The AE distribution for continuum and neutral generic B background 
are described by polynomials. The F distributions for the signal and the backgrounds are described 
by a Gaussian with different widths above and below the mean (a bifurcated Gaussian). 

The signal decay rate distribution of Equation 2 is modified to account for dilution coming from 
incorrectly assigning the flavor of Btag and is convoluted with a triple Gaussian resolution function, 
whose core width is about 1.1 ps. The decay rate distribution for B backgrounds is similar to 
that for signal. To account for their mis-reconstruction, the generic B backgrounds are assigned 
an effective lifetime instead of their respective measured B lifetimes. When evaluating systematic 
uncertainties, we allow for CP violation in the generic B background. This is described later in 
the text. The decay rate distribution for B^ J/ if) Kg is the same as that for signal and accounts 
for the known level of CP violation in that decay. The contimium background is modeled with a 
prompt lifetime component convoluted with a triple Gaussian resolution function. 

The results from the fit are 109 it 12 (stat) signal events, with S = —0.68 ± 0.30 (stat) and C = 
—0.21 lb 0.26 (stat). We also obtain for the aforementioned mutually exclusive tagging categories, 
the following numbers of continuum events: NLepton = 17 ± 5, NKaoni = 38 ± 8, NKaon2 = 101 ± 

12, NKaonPion = 102 ± 12, Npion = 115 ± 12, Nother = 94 ± 11 and N„oTag = 227 ± 17. 

Figure 2 shows the distributions of mESi and F for the data. In these plots the signal has 
been enhanced by cutting on |A^| < 0.1 GeV for the tties plot, jtzes > 5.275 GeV/c^ for the AE' 
plot and by applying both of these constraints for the F plot. After applying these requirements to 
the signal (background) samples that arc used in the fit, they arc reduced to a relative size of 83.1% 
(24.3%), 85.0% (21.1%) and 73.1% (2.8%) for the mEs, AE, and F distributions respectively. 

Figure 3 shows the At distribution for S° and B^ tagged events. The time-dependent decay- 
rate asymmetry [N{At) -T^{At)\/[N{At) + T^{At)\ is also shown, where N (TV) is the decay-rate 
for B^ (5°) tagged events and the decay-rate takes the form of Equation 2. 

4 Systematic Studies 

Table 1 summarises the systematic uncertainties on the signal yield, S and C . Each entry in the 
table indicates one systematic effect and these are added in quadrature to give the totals presented. 
These include the uncertainty due to the PDF parameterisation, evaluated by fixing both the signal 
and the background PDF parameters to their nominal values and varying them within uncertainties; 
the effect of SVT mis- alignment; the uncertainty due to knowledge of the Lorentz boost and z-scale 
of the tracking system, and knowledge of the event-by-event beam spot position. 

The uncertainty coming from the fit bias is estimated by performing ensembles of mock exper- 
iments using signal MC which is generated using the GEANT based BABMi MC simulation [15], 
embedded into MC samples of background generated from the PDF. The deviation from input 
values is added in quadrature to the error on the deviation in order to obtain the fit bias uncer- 
tainty. Most, but not all of the inclusive charmonium final states which dominate the generic B 
background, are precisely known from previous measurements. Their yields are then fixed in the 
fit. As a crosscheck, we allow the backgrounds to vary in the fit to data to validate the expected 
yields and to provide a systematic uncertainty. We also apply an additional systematic uncertainty 
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Figure 2: Signal enhanced distributions of niEs (top), AE (center) and F (bottom) for the data 
(black points). The (blue) solid line represents the total likelihood, the (red) dashed line is the 
sum of the backgrounds and the (black) dotted line is the signal. The undulations in the back- 
ground model are the result of limited MC statistics available for defining the two-dimensional 
non-parametric PDFs. 



to account for neglecting the small correlation between mss find AE in signal and neutral generic 
B background events. 

In order to evaluate the uncertainty coming from CP violation in the B background, we have 
allowed S and C to vary between -|-1 and —1 for the neutral generic B background, and for the 
direct CP asymmetry to vary between -1-0.5 and —0.5 for the charged generic B background. The 
CP parameters in B^ — > J ftp Ks are varied within current experimental knowledge [12]. 

The generic B background uses an effective lifetime in the nominal fit. We replace this with the 
B lifetime to evaluate the systematic error due to CP background lifetime. There is also a small 
asymmetry in the tagging efficiency between B^ and B^ tagged events, for which a systematic 
uncertainty is evaluated. We also study the possible interference between the suppressed b ucd 
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At(ps) 



Figure 3: The At distribution for a sample of events enriched in signal for (top) and (middle) 
tagged events. The dotted lines are the sum of backgrounds and the solid lines are the sum of signal 
and backgrounds. The time-dependent CP asymmetry (see text) is also shown (bottom), where 
the curve is the measured asymmetry. 



amplitude with the favored b — > cud amplitude for some tag-side B decays [16]. The difference in 
the distribution of F between data and MC is evaluated with a large sample of B ^ D* p decays. 

There are additional systematic uncertainties that contribute only to the branching fraction. 
These come from uncertainties in charged particle identification (5.2%), vr" meson reconstruction 
(3%), the J/ip branching fractions (2.4%), tracking efficiency (1.2%) and the number of B 

meson pairs (1.1%). The systematic error contribution from MC statistics is negligible. 

5 Summary 

The 109 lb 12 signal events correspond to a preliminary branching fraction of 
B{B^ J/V' vr°) = (1.94 ± 0.22(stat) ± 0.17(syst)) x 10"^ 
which is consistent with previous measurements from the B Factories. We determine the prelimi- 
nary CP asymmetry parameters 

C = -0.21 ± 0.26(stat) ± 0.09(syst), 

5 = -0.68 ± 0.30(stat) ± 0.04(syst), 
where the correlation between S and C is 8.3%. The value of S is consistent with SM expectations 
for a tree-dominated b — > ccd transition of S" = — sin 2/3 and C = 0. 
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Table 1: Contributions to the systematic errors on the signal yield, S and C. Additional systematic 
uncertainties that are applied only to the branching fraction are discussed in the text. 



Contribution 


Signal yield 


S 


c 


SVT mis-alignment 


— 


±0.002 


±0 


002 


Boost and z-scale 


+0.08 
-0.16 


±0.004 


±0 


001 


Beam spot position 




±0.007 


±0 


002 


PDF paramctcrisation 


+3.21 
-2.88 


±0.013 


+0 

-0 


009 
012 


Fit bias 


±3.00 


±0.030 


±0 


060 


Generic B background yields 


±3.52 


±0.003 


±0 


020 


Choice of ID v 2D PDFs 


±2.92 


±0.020 


±0 


002 


CP content of B background 


+0.40 
-0.26 


+0.020 
-0.018 


±0 


058 


CP background lifetime 


±0.67 


±0.010 


±0 


010 


Tagging efficiency asymmetry 


±0.02 


±0.000 


±0 


020 


Tag-side interference 




±0.004 


±0 


014 


Fisher data/MC comparison 


±0.70 


±0.004 


±0 


004 


Total 


+6.43 
-6.26 


+0.044 
-0.043 


±0 


093 
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